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I. INTRODUCTION 

Iron-based high- Tc superconductorai^i^i^i^i^i^i^ have a 
complex Fermi surface that is determined by several 
bands, an effect resulting from the hybridization of the 
3d orbitals of iron .^'^'^i-'^^'^^i^^ Band structure calcula- 
tions have shown that the bands that define the two 
hole pockets around the F point have mostly d^z and 
dyz character, while the two electron pockets around the 
M point have dxz, dyz, and a smaller amount of dj.y 
contributions j^i^'^i-'^^'^^i^^'^'^ For this reason it is important 
to understand superconductivity in multiorbital systems 
in general terms. 

Among the first to address this complex problem sev- 
eral years ago were Suhl et al^ using a model consist- 
ing of two orbitals, one s and one d, that did not hy- 
bridize with each other. Thus, in this case each band 
was determined by one single orbital. They showed 
that BCS pairinglS, could occur in each band and, since 
in the most general case the electron-phonon interac- 
tion would have different strengths for electrons in the 
different bands, it was proposed that two different su- 
perconducting gaps could arise. Almost 50 years were 
needed to observe experimental evidence of this phe- 
nomenon. In 2001, superconductivity with Tc = 39 K 
was observed in MgB2.— Despite the high- Tc, it became 
clear that the BCS mechanismlS was at play and for 
the first time two different superconducting gaps were 
0bserved-i2i2ai^i2^iS^i2ii2^ As shown in Ref. [li, the Fermi 
surface (FS) is determined by two bands: the tt band 
formed by the Pz orbitals of B, and the a band consti- 
tuted by a linear combination of the Px and Py B-orbitals. 
Although three orbitals determine the FS, it is interesting 
to notice that only two different BCS gaps are observed. 
This occurs because two of the three orbitals hybridize 
with each other and determine one single band, which 
couples strongly to the lattice phonons. This opens a 
large superconducting gap on the cr FS. The other or- 
bital, Pz, does not hybridize and forms the tt band that 
couples weakly to the lattice phonons determining a sec- 
ond, smaller, superconducting gap at the FS of the tt 
band. Thus, the number of different gaps that can arise 
in a multiorbital system is related to the degree of hy- 
bridization among the orbitals. Also note that in this 



early effort interband hopping of pairs of electrons be- 
longing to the same band was included but the possibility 
of interband pairing, i.e., pairs formed by electrons be- 
longing to two different bands, was not considered. 

In this paper, the subject of superconductivity in mul- 
tiorbital systems is revisited, in particular to shed light 
on the possible symmetry of the pairing operator of the 
pnictides superconductors. The motivation is that the 
pairing operators that have been discussed the most thus 
fari4^i27^i29^i3i^^^^^ assume that only intra- 
band pairing should occur, namely the two electrons 
of the Cooper pair belong to the same bandj^ How- 
ever, numerical simulations'^^ performed on a two-orbital 
mode l — 1^ '''^° for the pnictides favor an interorbital pair- 
ing operator that, when transformed to the band rep- 
resentation, results not only in intraband pairing but 
it includes interband pairing as well.**^ For the pnic- 
tides, interband hopping of pairs formed between elec- 
trons in the same band is often denoted as "Inter- 
band Superconductivity" 1^ The situation discussed in 
the present paper is different and involves Cooper pairs 
where the two electrons come from two different bands, 
which we will call "Interband pairing" . Using symmetry 
arguments and mean-field approximations, the plausibil- 
ity and physical meaning of such an interband pairing in 
multiorbital systems will be discussed. 

Interband pairing has previously been addressed in 
the context of Quantum Chromodynamics (QCD) and 
cold atomSf^i^ heavy fermions,*'^ cuprates*"', and BCS 
superconductivity.^^ In the case of heavy fermions, it was 
argued that interband pairing could occur if two Fermi 
surfaces arising from different bands are very close to 
one other)^ while in QCD and cold atoms it was pre- 
sented as a possibility for the case of sufficiently strong 
attractive pairing interactions, or for weaker attractions 
among particles with very different masses/^J^- As it will 
be discussed for a simple model in Sec. IIIIl three dif- 
ferent regimes, shown schematically in Fig. [l] can re- 
sult from a purely interband pairing as a function of the 
strength of the pairing potential g: (1) a normal regime 
where the ground state is not superconducting (namely in 
purely interband pairing an infinitesimal attraction does 
not lead to superconductivity), (2) an exotic supercon- 
ducting "breached" regime where gaps open at the nor- 
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mal Fermi surfaces while new Fermi surfaces defining re- 
gions containing unpaired electrons are created, and (3) a 
superconducting regime resembling BCS states, at large 
attractive coupling4^ 
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FIG. 1: Schematic representation of the three regimes that 
can arise as a function of the strength g of an interband 
pairing attraction. The label "normal" denotes a non- 
superconducting state. The case "breached" is an exotic 
regime with superconductivity and gaps, coexisting with 
Fermi surfaces (or several nodes) and electrons that do not 
pair. "BCS" is the large attraction region, where the ground 
state resembles that of a BCS superconductor and all elec- 
trons participate in the pairing. 

The paper is organized as follows. In Sec. [TTl the gen- 
eral form of pairing operators in multiorbital systems will 
be presented, remarking how the symmetry is determined 
by the spatial and the orbital characteristics of the op- 
erator. The interorbital pairing operator with B2g sym- 
metry obtained numerically in a two-orbital model for 
the pnictides is discussed, emphasizing that this oper- 
ator presents a mixture of intra and interband pairing 
in the band representation. In Sec. Illli a simple toy 
model with pure interband pairing attraction is intro- 
duced. This simplified model is discussed in order to 
illustrate the effects of interband pairing on observables, 
such as the occupation number and the spectral func- 
tions A(k, uj). The stability of the interband paired state 
is also discussed. The occupation number, spectral func- 
tions, and the stability of the i?2g pairing state are the 
subject of Sec. IIVI which is directly related to the physics 
of pnictides, while Sec. |V]is devoted to the conclusions. 

II. PAIRING OPERATORS IN MULTIORBITAL 
MODELS 

In single-orbital models, the symmetry of a spin-singlet 
pairing state is completely determined by the properties 
of its spatial form factor. More specifically, the pairing 
operator will have the form: 

A(k) = /(k)(ck,TC_k,i - Ck,iC-k,T), (1) 

where Ck.o- destroys an electron with momentum k and 
spin projection a and /(k) is the form factor that trans- 
forms according to one of the irreducible representations 
of the crystal's symmetry group. Thus, /(k) determines 
the symmetry of the operator. These form factors depend 
on the lattice geometry and generally they may be very 



complex. However, in materials with short pair coherence 
lengths, such as the high-Tc cuprates, the assumption 
that the two particles that form the pair can be very close 
to one other is usually made. The Cu-oxide planes in the 
cuprates have the symmetry properties of the D^h group 
and the case /(k) = cosfca, — cosfcy, which transforms 
according to the irreducible representation Big, provides 
the well-known d-wave symmetry pairing. 

In multiorbital systems, on the other hand, a spin sin- 
glet pairing operator will have both spatial and orbital 
degrees of freedom and it will be given by 

A(k) = /(k)TQ,^(dk,a,TC?-k,/3a - dk.Q.i'^-k./J.T), (2) 

where dk.a.o- destroys an electron with momentum k, in 
orbital a, and with spin projection cr, /(k) is the spatial 
form factor as indicated above, and Ta^fs is a matrix in the 
space spanned by the orbitals involved. The dimension 
of T is equal to the number of orbitals that are considered 
to be of relevance. In this case, notice that the symmetry 
of the pairing operator would in general be determined 
by the product of the symmetry properties of /(k) and 
the symmetry of the orbital contribution Ta.p. Only if 
Tq./3 is the identity matrix do the orbital contribution 
become trivial, because the identity matrix transforms 
according to Aig. Thus, this is the only case where /(k) 
fully determines the symmetry of the pairing operator, 
as in the single-orbital example. 

The minimum model for the pnictides considers 
the two orbitals dxz and dyz, which are strongly 
hybridized. '^^'■^^ All the possible pairing operators, up 
to nearest-neighbor distance, that are allowed by 
the lattice and orbital symmetries have been al- 
ready calculated i'^''''^^i"^^'^°i^^'^^ Numerical simulations 
performed on the two-orbital model suggest that the fa- 
vored pairing operator at intermediate couplings, where 
the state is both magnetic and metallic,'^®'''*' has symme- 
try B2g and is given by Eq. ^ with /(k) = (cosfca, -|- 
cosfcy), which transforms according to Aig, and t — ai 
which transforms according to B2^ (where ai are Pauli 
matrices). Thus, the non-trivial symmetry under rota- 
tions arises from the orbital portion of the operator. This 
pairin g op erator has been studied at the mean-field level 
in Ref. ®. In the orbital representation, the Bogoliubov- 
de Gennes Hamiltonian matrix is given by: 
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with 

i,xx = — 2t2 cos kx — 2ti cos ky — At^ cos k^ cos ky — 

iyy ~ —2ti cos kx — 2t2 cos ky ^ Ats cos kx cos ky — ^, 

ixy = -UiSirikxSmky, (4) 

and 

Ak = V^(cos kx + cos ky), (5) 
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where V = VqA, with Vq being the strength of the pair- 
ing interaction and A the mean-field parameter obtained 
by minimizing the energy. Since the two orbitals are hy- 
bridized via S,xy, in the band representation the Hamil- 
tonian matrix becomes: 



— 



ei Vb Va 

62 -Va Vb 

Vb -Va -£2 

Va Vb -ei 



(6) 



where Va and Vb are given by 

Va = 2?/(k)w(k)Ak, 

Vb = iv(k)' - u{k)')A^, 



(7) 



(8) 



and M(k) and u(k) are the elements of the change of basis 
matrix U given by 



U 



/u(k) v{k) 

v{k) -u{k) 

v{k) u{k) 

V -u{k) v{k) 



(9) 



with U^^ = . Remember that Va and Vb, are func- 
tions of the momentum k, and u(k)^ + u(k)^ = 1. Thus, 
it is clear that in the band representation, in addition 
to the intraband pairing given by Va, there is also in- 
terband pairing given by Vb- Among pairing operators 
compatible with the symmetry of the model, the ones 
that do not lead to interband pairing not only cannot 
mix orbitals, but have to contain t — ao, i.e., the iden- 
tity matrix-^*' One example for such an operator would 
be the s± pairing 

The discussion above describes general properties of 
hybridized multi-orbital systems. If the orbitals are hy- 
bridized, but not related to one another by symmetry, 
there is no reason to expect that the coupling between 
the electrons in each orbital, and the interaction that pro- 
duces the pairing, will have the same strength for all the 
orbitals and lead to a unit matrix in the orbital sector. 
Then, it is expected that interband pairing will arise in 
general and, thus, it is important to understand its con- 
sequences, providing the main motivation for the present 
manuscript. 



III. INTERBAND PAIRING 



A. Generic properties 



Model and non-interacting limit 



To address qualitatively the issue of interband pairing, 
postponing the matters of stability to the next subsec- 
tion (Sec. IIIIB]) . let us consider the following two-bands 



simplified model with interband pairing: 

a,(7 oi^fi 

(10) 

where a, (3,— 1,2 label two bands that are not hy- 
bridized, a is the spin projection, and for simplicity 



ea(k) 



2ma 



ill) 



which gives parabolic bands that are degenerate at k = 
with energy C, and with a chemical potential /i = 0. 
This can be considered as a crude representation of the 
two hole-pocket bands around the T point in the pnic- 
tides, but more importantly presents a simple toy model 
where the effects of interband pairing can be studied. As 
before, the parameter V = VqA is the product of an at- 
tractive potential Vq between electrons in the two differ- 
ent bands and a mean-field parameter A determined by 
minimizing the total energy. The band dispersion with- 
out the interaction is presented in the inset of Fig. [21 
The Bogoliubov-de Gennes matrix expressed in the basis 
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FIG. 2: (Color online) Mean-field band dispersion for the 
model defined by Eq. (|10p . for the indicated values of V 
(defined in the text) as a function of the momentum k — 
A^/ fc§ + ky. The case shown is for mi — 1, m2 — 2, and 
C — 2. Inset: Non-interacting band dispersion for the same 
parameters. 



expanded by B12 = {cj,.^ |, c_k,2,i; Ck_2,|j c_k,i,x} has the 
form: 




H 



This matrix can be diagonalized becoming 



Hr 



(12) 
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-Ea 



(13) 



4 



in the basis expanded by 

{7k,A.T ' 7-k,Ba , 7k3,T ' T-k.A,i} and 
basis matrix is given by 



U = 



Uk 

— Uk Uk 

Mk t'k 

— Uk Itk 



Bab 
the change-of- 



(14) 



with ut 



The four energy eigenvalues are 



S,(k) = ± 



£2 



)2 + l/2 



(15) 



where the first positive (negative) sign corresponds to 
the eigenvalues of the upper (lower) block, labeled Ea 
and -Eb [Eb and -Ea) in Eq. ^ and in Fig.H The 
second sign differentiates between the two solutions in 
each block. Figure [2] shows the eigenvalues for V = and 
V = 0.5. When V = 0, then Ea = £2 and Eb = ci, the 
two bands define two circular Fermi surfaces with radius 
kpi and kF2 where they cross the chemical potential {p = 
0). This is illustrated schematically in Fig.[3l 



unpaired 
electrons 1 




paired 
electrons 



while in basis Bab the number operator is (/=A,B) 

nz(k) = V7i,,7k,i,.. (17) 



The total electronic occupation of the system is given by 
n(k)=^n„(k)=^nKk). (18) 

a I 

Then, for V ^ wc find that n{k) = for |k| < |ki.i|, 
ri(k) = 2 for |ki?i| < |k| < \kp2\ since in this region 
ni(k) = n^(k) = 2 while n2(k) = nB(k) = 0, and finally 
n(k) = 4 for |k| > \\iF2\ since both orbitals are totally 
filled with electrons. These results are represented by the 
dashed lines in Fig. [H 
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FIG. 4: (Color online) Mean-field state population as a func- 
tion of momentum along the diagonal = ky for (a) the 
whole system; (b) band B, (c) band A; (d) orbital 1; and (e) 
orbital 2 for the indicated values of the pairing potential V, 
and mi = 1, m2 = 2, and C — 2. 



2. Weak attraction 



FIG. 3: (Color online) Schematic diagram of the FS deter- 
mined by the two parabolic bands of the simple model used 
in Sec. IIIII kpi and kF2 indicate the Fermi momentum of the 
two bands when V = 0, while k'p^ and k'p2 indicate the posi- 
tion of the Fermi momenta for the case of a finite but small 
pairing potential. The shaded rings indicate the regions with 
width Ai and A2 in momentum space where electrons can 
pair. The white region in between the rings contains unpaired 
electrons in band 1. 

Notice that number operators can be defined in the two 
bases that are being considered here, i.e. B12 and Bab, 
which of course become equivalent when V=0. Thus, in 
the basis B12 the number operator is, 

"a(k) = 4,a,<TCk,a,a, (16) 

a 



In the nontrivial case of V different from zero, the 
bands ±Ea and ±_Eb result from the hybridization of 
the ei and £2 bands due to V. It is interesting to ob- 
serve that an internal gap opens at the crossing of bands 
Ea with —Eb above the chemical potential and between 
Eb and —Ea below the chemical potential, as indicated 
with circles in Fig. [5] where results for V = 0.5 are dis- 
played. These are very important differences with re- 
spect to conventional BCS calculations were all the action 
is restricted to the original Fermi surfaces: for multior- 
bital models, gaps can open in other portions of the band 
structure as well. 

While the ±Ea bands are separated by a gap, the 
bands ±Eb still cross the chemical potential determin- 
ing two Fermi surfaces at k'p^ > kpi and at k'p2 < kp2, 
even in this pairing state (again, this is different from the 
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one-orbital pairing standard BCS ideas). The new Fermi 
surfaces are shown schematically in Fig. [31 the interior 
FS has expanded and the exterior one has contracted. In 
fact, it will be shown below that the effect of V is to try 
to equalize the two Fermi surfaces, as this pairing attrac- 
tion grows in magnitude. Calculating 7i(k) we observe 
that nCk) = for k < k'p^ and k > k'p2, where is the 
element of the change-of-basis matrix in Eq. This 
agrees with the BCS expression for n(k) but it jumps 
discontinuously to n(k) = 2 for /c^^ < k < k'p^- Such a 
result is shown by the solid lines in Fig. [31^ a). The jumps 
indicate the existence of the two Fermi surfaces, which 
are here present even in the paired state. Thus, some 
electrons in the region in between the two Fermi surfaces 
may behave like normal unpaired electrons. 

A better understanding of the electronic behavior can 
be achieved by studying the electronic density in the two 
bases Bab and B12. We find that n^(k) = 2w^, as shown 
in Fig.m^c), which is the standard BCS behavior in agree- 
ment with the fact that bands ±Ea are separated by a 
gap. Thus, all the electrons in this band participate in the 
pairing and they do not have a FS. On the other hand, it 
can be shown that ^^(k) = ?T.A(k) = for k < k'p^ and 
k > kp2, but for k'p^ < k < k'p2 there is a discontinuous 
change of behavior to nB(k) = = I — [Fig. Hllb)]. 
Thus, the Fermi surfaces are determined by electrons in 
the band Eb- However, note that in this region nyi(k) is 
an increasing function of |k| while ^^(k) is a decreasing 
function of |k| which satisfies n^(k) + riB(k) = 2 for all 
|k|. 

This behavior can be better understood by calculat- 
ing the (photoemission) spectral functions j4(k, w), which 
allow us to obtain n(k) = A{'k,ijj)duj. In the non- 
interacting case, shown in Fig.O^a), the spectral function 
shows two peaks, corresponding to the two bands ei and 
£2, for each value of |k|. This is the expectation for free 
electrons in a non-interacting multiorbital system. The 
two non-interacting FS's are located where each band 
passes across the chemical potential. 

When the pairing interaction becomes finite, the spec- 
tral function develops four peaks for each value of |k| as 
it can be observed in Fig. Efb) for V = 0.5. This is also 
the expected result since the BCS interaction generates 
a "shadow" or Bogoliubov band for each band present in 
the non-interacting system. For example, close to |k| = 
the bands Ea and Eb have almost all the spectral weight, 
given by which is very close to 1 in this region, and fol- 
low a dispersion similar to the non-interacting bands ei 
and €2, while the bands —Ea and —Eb appear with very 
small spectral weight given by = 1 — m^. The latter 
are the Bogoliubov or "shadow" bands. These shadow 
bands appear above the chemical potential for large val- 
ues of |k|, as expected. 

But what happens in the intermediate region k'p^ < 
k < fc^2? We see that the band Eb crosses with —Eb 
at k'p^ determining a FS. Most of the spectral weight 
below the chemical potential belongs to the band —Eb 
which contains the unpaired electrons but the band —Ea 
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FIG. 5: (Color online) (a) Spectral functions A(k,uj) for the 
non-interacting case, V = 0. (b) Same as (a) but for V = 0.5. 
Note that there are still Fermi surfaces in this intermediate 
coupling case, namely gaps open far from the chemical po- 
tential but not at the Fermi level, (c) Same as (a) but for 
V — 1.0. The presence of a gap for all the momenta shown is 
now clear. 



has appreciable shadow spectral weight indicating that 
there are also some paired electronic population as in- 
dicated in Figs. Sfb) and (c). As k increases, spectral 
weight is transferred continuously from —Eb to —Ea, 
behavior associated with the internal gap opened by the 
pairing interaction, so that when k approaches k'p2 most 
of the spectral weight below the chemical potential is in 
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—Ea, paired electrons, and —Eb, unpaired electrons, has 
shadow spectral weight as seen in Figs-Sfb) and (c). At 
k — the second FS is determined by the crossing of 
±Eb at the chemical potential indicated by the sudden 
jump in n^. Thus, the unpaired electrons in this region 
coexist with paired electrons and the spectral functions 
do not resemble the non-interacting ones. 

It is also illuminating to analyze what happens with 
7T.(k) in the basis Bi2. While ni(k) — 7i2(k) = 2v^ in the 
regions where there are no unpaired electrons, we find 
discontinuities associated with Fermi surfaces in the two 
distributions that are given by 712 (k) = and ni (k) = 2 
for kpi < k < k'p2 [Figs. IHd,e)]. This indicates that 
the pairing interaction V has been able to promote some 
electrons from above to below kp2 in orbital 2. Also, 
electrons have been transferred from their original loca- 
tion in orbital 1, in the neighborhood of kp2 and above 
kpi, to both orbitals 1 and 2 around kpi. These are 
the electrons in 1 and 2 that have become paired (the 
pairing is indicated by the shadowed circular regions in 
Fig. [3]). But the interaction was not strong enough to 
provide pairing partners to all the extra electrons origi- 
nally in orbital 1 and, thus, they have been left unpaired 
in between the two paired regions, as indicated in Fig.O 

Thus, the interband pairing attraction creates pairs of 
electrons belonging to different orbitals within an interval 
l^kpi around each of the two original Fermi surfaces. The 
width of the pairing region increases with V . The pair- 
ing partners are obtained by promoting electrons with 
momentum k k, kp2 and k w kpi in both orbitals and 
by moving electrons from the more populated to the less 
populated orbital. This creates the conditions to pair 
electrons near both Fermi surfaces. The electrons in band 
1 that could not find promoted partners remain unpaired. 
Whether this state is stable or not depends, of course, on 
the balance between kinetic and pairing energies, which 
will be discussed in Sec. IIII Bl 



3. Strong attraction 
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FIG. 6: (Color online) Mean-field state population as a func- 
tion of the momentum (main diagonal) for (a) the whole sys- 
tem, (b) band 1, and (c) band 2 for the indicated values of 
the pairing potential V , and for mi = 1, m2 = 2, and C = 2. 
The case V = \S) illustrates the "strong" pairing attraction 
regime where gaps open in the original Fermi surface. 



Thus, notice that while band 1 contained more elec- 
trons than band 2 for V — Q [see Figs. [Bl^b) and (c)], 
the interband pairing mechanism transfers electrons from 
one band to the other so that both bands have the same 
number of electrons in the superconducting state. Con- 
sequently, the smaller FS expands and the larger one 
shrinks. Then, when the pairing becomes strong enough, 
the two Fermi surfaces become equalized and no unpaired 
electrons remain. Whether this situation can be achieved 
will depend on the strength of the interaction and the en- 
ergy balance, as discussed in the next subsection. 

If the two non-interacting Fermi surfaces are very close 
to each other in momentum space, even relatively weak 
pairing interaction could effectively be strong enough to 
make the interband pairing resemble BCS pairing as in 
the case of large V in the present example. 



As the interaction V increases further, the number 
of unpaired electrons is reduced. This means that k'p-y 
and k'p2 become closer to each other making the size of 
the intermediate region with unpaired electrons in Fig. [3] 
smaller. Eventually, the two momenta become the same 
k'p^ = k'p2 for V = 0.71, and for V > 0.71 the region 
with unpaired electrons vanishes and a full gap opens in 
the system whose physics now resembles BCS, except for 
the fact that the pairs are constituted by electrons from 
different orbitals. The electronic population of the sys- 
tem in such a case, e.g. at ^ = 1.0, is presented in Fig. [6] 
and the corresponding spectral functions are shown in 
Fig. ^c). It can be shown that now ni(k) = n2(k) for 
all values of k and all the particles around the two non- 
interacting Fermi surfaces now participate in the pairing. 
The spectral functions show four peaks, i.e. Bogoliubov 
bands for all values of k, as it can be observed in Fig.[5l^c). 



B. Stability of the interband paired state 

1. The case without intraband pairing 

As discussed in the Introduction, the possibility of in- 
terband pairing has been previously discussed in the con- 
text of QCD and cold atomic matter. Similar effects on 
the FS as found in the present study were described, al- 
though the physics was different because in the QCD 
context each band contained different kinds of particles 
and the pairing was thus not able to promote particles 
from the majority to the minority band^ii^ The issue 
of stability was explored in the QCD framework, and it 
was found that a purely interband paired state could be 
stabilized for pairing attractions above a certain cut-off 
value, which could become very small for a large differ- 



ence between the masses of the two paired species i^^i 

In the case of our model, however, we have found (see 
below) that the purely interband-paired state only be- 
comes stable when the attraction is sufficiently strong 
that no unpaired particles are left, i.e. when the two 
shaded regions overlap and the unpaired region in Fig. [3] 
vanishes. This means that, although the pairs would 
be formed by electrons in different orbitals, the physics 
would be analogous to BCS. A gap will be opened in the 
full Fermi surface of the simple model studied here. 

In order to study the issue of stability, let us assume 
that the interaction term responsible for the intcrband 
attraction is given by 

k,k',Q: 

(19) 

where Vk,k' = —Vq and A'' is the number of sites. Per- 
forming the standard mean-field approximation: 6k' = 

(c-k',-Q,iCk',Q,T) and = (Ck.Q,T'^-k.-Q,i^ ^^'^ making 
the substitution c[_„^^cli^__^^^ = bl + (4,„jcLk,_„^x " 

6jj.) (and an analogous substitution for the product of an- 
nihilation operators), the mean- field results are obtained. 
As usual, the fluctuations around the average given by 
("^k a T'^-k -a I ~ ^k) ^re assumcd to be small. Defining 
A = -i- bk — j^ J2k ^k obtain the following mean- 
field Hamiltonian: 

HmF = ea(k)c]^^^Ck,a,o- 

a,cr 

-V,A ^ (c^^^cLk,^,^ + h.c.) + 2V,A'N. (20) 

Equation (|12p can be recovered by defining — VqA = V, 
and disregarding the constant last term of Eq. (|20p . We 
can calculate the total energy Emf for Eq. (|20p as a func- 
tion of = VbA. If for a given V the energy has a 
minimum, this indicates that the interband-paired state 
is stable. Note that having a term linear in A in Eq. ((20|) 
is not sufficient to conclude the appearance of a supercon- 
ducting state at small A, since the sign of the coefficient 
of the linear term can change sign with V. A similar 
situation occurs for the magnetic state of undoped pnic- 
tides: a finite Hubbard U must be reached to stabilize 
the "striped" stated. Rctur ning to superconductivity, 
there are two regions of interest: (i) < V < 0.71, which 
corresponds to the case in which two Fermi surfaces are 
present in the paired state; (ii) V > 0.71, which corre- 
sponds to the fully gapped case. In Fig.[7Ka), Emf/N vs. 

V is shown for different values of Vq . It can be observed 
that a second minimum develops for Vq > 2 and it be- 
comes stable for Vq > 3. The minimum always occurs for 

V > 0.71 which means that it corresponds to the case in 
which there are no unpaired electrons in the system. The 
results shown in the figure are robust in the sense that 
changes in the values of mi, m2 or the chemical poten- 
tial were not found to stabilize the state with unpaired 
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FIG. 7: (Color online) (a) Mean-field energy Eq. (f20|) per site 
vs. V = VJ)A for different values of Vq. The case V = 0.71 
is indicated with dashed lines: it separates the regions corre- 
sponding to nodal and nodeless states, (b) Mean-field energy 
(Eq. I20|l per site with intraband pairing with strength /2 
vs. 1/ = VoA, for different values of Vb. 

electrons. Thus, in this respect an attraction that is only 
interband can only lead to a stable superconducting state 
in the strong attraction region. 



2. Stability when both inter and intraband pairing coexist 

The results of the previous paragraphs may seem neg- 
ative with respect to the relevance of the "intermediate" 
state with simultaneous coexistence of pairing and Fermi 
surfaces. However, as pointed out in Sec. [TTl most of the 
pairing operators allowed by the lattice and orbital sym- 
metry in the pnictides are characterized by a mixture of 
both intra and interband pairing. Thus, it is important 
to consider such a situation in our simple model as well. 
In the case of the i?2g pairing operator, Eqs. ([6]), (O, and 
([5|) indicate that the pairing is purely interband only for 
kx = or TT and fcy = or tt, because Va — along 
these lines. Thus, let us now consider our simple model 
in a Brillouin zone (BZ) defined by — tt < kx,ky,< tt 
and with the addition of intraband pairing with intensity 
Vo/2 that is vanishing for fc^; = or tt and ky = ov 
TT, i.e., the pairing is purely interband only along those 
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directions. The energy bands now behave as in Fig. [2] 
only along (0, 0) — (tt, 0) and (0, 0) — (0, tt), while the dis- 
persion along any other direction is shown in Fig. [8] for 
y = and y = 0.5. 



(7t,0) 




FIG. 8: (Color online) Mean-field band dispersion for the 
model given by Eq. (|10|) (along the main diagonal) for the 
indicated values of the pairing potential V, with the addition 
of an intraband pairing with strength Vo/2, as described in 
the text. 




Performing the mean-field approximation similarly as 
explained above, we have found that the superconduct- 
ing state now becomes stable on both sides of the original 
critical value V = 0.71. Figure [7]Jb) shows that the pair- 
ing state is stabilized for Vb > 2.5 but the value of A 
where the minimum is located is such that V < 0.71 
and, thus, two nodes will be present along the x and y 
axes. Increasing the value of Vb, the minimum eventually 
occurs for V > 0.71. For these larger values of Vb, there 
would consequently be no nodes. 

Then, in this section it has been shown using a simple 
model that the interorbital paired state can become sta- 
ble if the attraction Vb is sufficiently strong. In this case, 
it is the nodeless case that is stable even for purely inter- 
band attraction. In addition, the very interesting novel 
phase with coexisting nodes and unpaired electrons in 
the majority band also requires intraband pairing to be 
stable, with strength similar to that of the interband, at 
least in parts of the BZ. 



IV. THE INTERORBITAL Bag PAIRING 
OPERATOR 

In the previous section, a simple model was presented, 
both with exclusively interband pairing and with both 
inter and intraband pairing, and it was found that the 
intraband pairing stabilizes the state with a mixture of 
superconductivity and metallicity. As mentioned in the 
Introduction, it is expected that in the most general cases 



FIG. 9: (Color online) (a) One-particle spectral function for 
the two-orbital model Eq. Q, with vanishing pairing interac- 
tion V — 0. Parameters: ti — 1.3, t2 — —1, ts — ti = —0.85, 
fj, = 1.54. (b) Same as (a) but for V = 0.5. 



the pairing operators allowed by the symmetry of the lat- 
tice and of the orbitals will, in the band representation, 
have both intra and interorbital pairing. Thus, now we 
will present and discuss, at the mean field level, the oc- 
cupation number and the spectral functions for the pair- 
ing operator obtained from the numerical study of the 
two-orbital model for the superconducting state of the 
pnictidcs introduced in Sec. |TT1 



A. Non-interacting limit 

Let us start with the non-interacting case in which 
V^ = 0. The spectral functions along high-symmetry di- 
rections in momentum space are presented in Fig. 9(a) 
and, as expected, they reproduce the non-interacting 
band dispersion. We also show the total occupation 
number n(k) along the same directions [dashed lines in 



Fig.fTUTa)]. as well as the occupation number for each or- 
bital nx{k) [dashed hues in panel (b)] and ny{k) [dashed 
lines in panel (c)] and the orbital occupation in the 
quadrant of the first BZ defined by < fc^:, fey < tt in 
Fig. [TIT a) . It is clear that the electron-like and hole-like 
Fermi surfaces are determined by an admixture of the 
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(a) 




(3t,it) (0,0) (jt,0) 



FIG. 10: (Color online) (a) Total occupation number n(k) for 
the two-orbital model Eq. (|3]) with pairing interaction V = 0.5 
(continuous lines) and V = (dashed lines). Parameters: 
ti = 1.3, t2 = —1, ts = t4 = —0.85, jj. = 1.54. (b) Same as 
(a) but for the orbital d^z- (c) Same as (a) but for the orbital 
dyz. (d) Same as (a) but for band 1. (e) Same as (a) but for 
band 2. 



two orbitals. 

On the other hand, in the band representation, band 
1 determines the electron pockets while band 2 forms 
the hole pockets as it can be seen from the behavior of 
ni(k) and n2(k) (indicated by the dashed Unes in panels 
(d) and (e) of Fig. [TO]) and by the light (orange) and 
dark (red) surfaces in Fig. [TlTb) where the FS is also 
indicated. It is clear that the electronic occupation of 
band 1 is smaller than the electronic population of band 
2, so that unpaired electrons would be expected to belong 
predominantly to band 1, as in the simple model of the 
previous Sec. IIIII It is interesting to notice that in the 
orbital representation, on the other hand, the electrons 
are equally distributed among the xz and yz Sd orbitals. 



B. Nonzero pairing 

Let us discuss what occurs when the pairing interac- 
tion becomes nonzero. To simplify the discussion, define 
the following points in momentum space: X = (tt, 0), 
Y = (0^, r = (0,0), and M = (7r,7r). As remarked 
in Ref. |40, the B2g pairing operator always has nodes 
along the X — Y direction because the spatial form fac- 
tor /(k) — coskx + cos ky vanishes along that line. But, 
as soon as V is finite, a gap opens along the T — M direc- 
tion [notice that along this direction the pairing is purely 
intraband since in Eq. ([S]), = and thus Vb = in 





0,5 

k /i 

y 







' 0,5 



FIG. 11: (Color online) (a) Occupation number n(k) for the 
two-orbital model Eq. ^ without pairing interaction for the 
orbital dxz (orange/light) and the orbital dyz (red/dark). The 
'floor' indicates the FS in red. (b) Same as (a) but for band 
1 (orange/light) and band 2 (red/dark) 



Eq. dS])]. Along T-X,r-Y, X~ M, and Y-M nodes 
associated to the different number of electrons in band 1 
and band 2 remain [notice that along these directions the 
pairing is purely interband since Va — because Eq. ([3]) 
and Eq. ^ become identical to one other]. When the 
pairing interaction V becomes strong enough to make 
ni(k) = 77.2 (k), as described in the simplified model pre- 
sented in the Sec. IIII A 31 these nodes vanish. Along any 
other direction in the BZ a mixture of intra and interor- 
bital pairing will be present.— 

Let us first consider a relatively small pairing V = 0.5. 
In Fig. [TOFa). n(k) is presented along high symmetry di- 
rections. Along Y — X there is no pairing and, thus, ?T.(k) 
is unchanged from the non-interacting case shown in the 
figure with dashed lines. Along X — M , where only inter- 
band pairing occurs, no effects are observed at the elec- 
tron pocket FS but a rounding in n(k) indicating pairing 
is observed at the hole pocket FS. However, n(k) shows 
discontinuities at two points indicating the existence of 
nodes. Along the diagonal direction M — T, where all 
the pairing is intraband, it is found that n(k) exhibits 
standard BCS behavior at both hole Fermi surfaces indi- 
cating the opening of gaps. Finally, along T — X it can be 
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observed a rounding of n(k) at the hole Fermi surfaces, 
indicating pairing, and a sharp jump at the electron FS. 

We can further analyze the pairing in the orbital rep- 
resentation. The occupation number for the orbitals xz 
and yz is shown in Figs. [TUTb-c). Along Y — X, where 
the pairing is zero, we observe how the FS for the elec- 
tron pocket at y {X) is totally determined by electrons 
in the orbital xz {yz) and how the population of each 
orbital varies smoothly between the two Fermi surfaces, 
always satisfying nxi}i) + ny{k) = 2. From X to M, 
nx(k) — ?ij^(k) ior kph < k < M indicating that the elec- 
trons at the hole pocket FS are paired in a wide region 
around it, but the pairing region is very narrow around 
the electron pocket because the pairing is reduced by 
the small value of /(k). Along the diagonal, i.e. from 
M to r, standard intraband pairing occurs at both hole 
Fermi surfaces and, thus, ria;(k) — ny{k) while from F 
to X n2;(k) = nj,(k) for T < k < kph indicating pairing 
around the hole FS and almost no pairing occurs at the 
electron pocket FS. The behavior n2:(k) — nj,(k) = 
for all k for which pairing occurs and ^^(k) unchanged 
from the non-interacting value for all other k, was ob- 
served for all values of V . For this reason, figures for 
n.i(k) in the orbital representation will not be shown for 
the additional values of V discussed below. 

The population of the different bands is presented in 
FigHni (d-e). The Fermi surfaces along Y — X are clearly 
determined only by band 1, which is the band that forms 
the electron pockets. It can also be observed that there is 
almost negligible pairing at the electron FS along X — M, 
but there is clear interband pairing at the hole FS. This 
is an indication that, due to the spatial variation of the 
pairing interaction, the attraction is much stronger at the 
hole pockets than at the electrons pockets. Also notice 
that at the hole pocket FS the pairing is interband and 
thus ni(k) — n2(k), but this does not happen along the 
diagonal direction M — T where intraband pairing occurs, 
and there are more paired electrons belonging to band 2 
than to band 1. Along the direction X — T again we 
observed a stronger pairing effect at the hole FS than at 
the electron one. 



and one finds a Fermi surface on both the hole and elec- 
tron pockets, indicating that V = 0.5 corresponds to the 
"breached" phase in Fig.[TJ However, the spectral weight 
determining the hole-pocket node is weak and might thus 
be missed in the analysis of experiments. The node on 
the electron pocket FS, on the other hand, is robust and 
it should be observed, if present. The same occurs along 
X — M; the signal for the node at the electron pocket FS 
should be robust while the one at the hole pocket FS will 
be weak. 



S. Intermediate attractive coupling 

As the pairing interaction increases, the nodes result- 
ing from the interband pairing should get closer to each 
other, as discussed in Sec. IIIII Figures [121 a-c) show 
n(k) ior V = 3 along high-symmetry directions in the 
band representation. One can see that there are still un- 
paired electrons, and V = 3 consequently falls into the 
"breached" region schematically represented in Fig. [TJ 

It is interesting to notice that while k'pf^ > kph, on the 
other hand k'p^ ~ kp^ indicating that the reconstruction 
around the hole pockets is much larger than around the 
electron pockets. This can also be observed in Fig.fTWb). 
where we observe unpaired electrons along X — M, but 
not along T — X. The effect occurs in part due to the 
smaller value of /(k) at the electron pockets but also be- 
cause, due to the band dispersions, the price in kinetic 
energy for interband pairing is much larger at the electron 
pockets than at the hole pockets. The spectral density 



in Fig. 13(a) further illustrates that the pairing interac- 
tion is more effective along T — X than along X — M. 
Only shadow spectral weight crosses the chemical poten- 
tial along T — X, while strong spectral weight crosses the 
chemical potential twice along X — M and leads to two 
nodes. 



3. Strong attractive coupling 



C. Spectral functions j4(k, u) 

In this section, we discuss the form of the spectral func- 
tions A(\i,Lu), which can be measured in angle resolved 
photoemission spectroscopy (ARPES) experiments, for 
weak to strong interorbital pairing. 



1. Weak attractive coupling 

The spectral function for V = 0.5 is depicted in 
Fig. |9(b)[ As discussed above, the interorbital B2g oper- 
ator leads to mtra-band coupling along the F — M line, 
where one consequently clearly sees the hole pockets to 
be gapped. Along T—X, the pairing is purely mter-band. 



Finally, let us consider a much stronger value of the 
pairing, such as ^ = 6, in the BCS region of Fig. [U for 
which the only nodes observed are those along X—Y, due 
to the vanishing of the pairing operator. Figures [T^ d-f). 
show that n(k) is discontinuous only along X — Y, while 
it is smooth along all the other directions indicating pair- 
ing. Note that ni(k) = n2(k) along T - X and X - M 
where interband pairing occurs, while along T — M, ni(k) 
is smooth but different for each band because the pair- 
ing is intraband. The behavior of the spectral functions 
displayed in Fig. 13(b) shows that spectral weight only 
crosses the chemical potential along the X — Y direction. 
In this situation, in the folded BZ, nodes should occur 
only at the points where the two electron pockets cross 
with each other, as indicated in Ref. 40. In the rest of 
the BZ an anisotropic gap will be observed. 
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FIG. 12; (Color online) (a) Total occupation number n(k) for 
the two-orbital model Eq. ((3} with pairing interaction V = 3 
(continuous lines) and V = (dashed lines). Parameters: 
ti = 1.3, t2 — —1, ts = t4 — —0.85, n — 1.54. (b) Same as 
(a) but for band 1. (c) Same as (b) but for band 2. (d) Total 
occupation number n(k) for the two-orbital model Eq. ^ 
with pairing interaction V — 6 (continuous lines) and V — 
(dashed lines), (e) Same as (d) but for band 1. (f) Same as 
(e) but for band 2. 
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FIG. 13: (Color online) (a)Spectral density for the two- 
orbital model Eq. (|3J for a pairing strength V = 3. The 
rest of the parameters are as in Fig. [O] (b) same as (a) but 
for strong pairing V — 6. 



form 

Vk,k' = V*(cos kx + cos fcj^)(cos k'^ + cos fc^), 
and that 

At(k) = A(k) = V*A{cosk^+cosky), 



(21) 



(22) 



where 



A(k) = -^Vu,k'(&k') 

k' 

At(k) = -5]yk,k'(&L) 



(23) 



D. Stability of the B2g pairing state 

Finally, let us discuss the important issue of the stabil- 
ity of the B2g pairing state. It will be assumed, follow- 
ing the notation in Appendix A of Ref. [i^, that Eq. ^ 
has arisen from an interorbital attractive potential of the 



V in Eq. (O is then given hy V = V* A and the mean- 
field energy Emf can be calculated for a given V* as a 
fimction of V. The results are presented in Fig. [Ml In 
this figure, a minimum for V ^ can be seen for V* > 3. 
This indicates that if the pairing attraction overcomes 
a finite threshold level, the mean field results presented 
here with V > 0.2 will be stable. Thus, if the attraction 
exists, both multinodal (breached) and states with nodes 
only along X-Y (strong coupling) are possible. 
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FIG. 14: (Color online) Mean field energy per unit site for the 
interorbital pairing with symmetry B2g for different values of 
the attraction V* as a function of = V*A. 

V. CONCLUSIONS 

Summarizing, in this manuscript the possibility of 
intra and interband pairing in multiorbital systems 
has been discussed. While interband pairing has pre- 
viously been studied in the context of QCD, cold 
atomsj^ii^ heavy fermionS)^ cupratea^, and BCS 
superconductivity;^ most of the pairing operators pro- 
posed for the pnictides are based on the premise that 
the pairing has to be purely intraband, i.e., both elec- 
trons in the Cooper pairs belonging to the same band, 
compatible with the assumption that in these supercon- 
ductors all the action must occur at the Fermi surfaces. 
However, symmetry considerations show that in models 
for the pnictides interorbital pairing is allowed with the 
two members of the Cooper pair belonging to different 
bandsji^ This is not surprising since the bands that de- 
termine the electron and hole Fermi surfaces consist of 
hybridized orbitals. These interorbital pairing operators 
give rise to not only intra but also interband pairing when 
the band representation is used. In addition, numerical 
calculations in a minimal two orbital model for the pnic- 
tides favor one of these non-trivial pairing operatorsi^^i^ 
As a consequence, a clear discussion of the role of inter- 
band pairing is necessary. 

The explicit calculations shown here of the electronic 
occupation in the orbital and the band representations, 
as well as the calculation of the spectral functions show- 
ing the distribution of Bogoliubov bands, may offer guid- 
ance in the interpretation of ARPES experiments. In 
particular, most ARPES measurements determine the 
FS in the normal state and study the opening of the su- 
perconducting gap by monitoring A(kp, lu) as they lower 
the temperaturci^i^i^iSi^ Notice that this approach 
would miss the nodes associated with the "breached" 
phase since in the superconducting state a gap would be 
observed in A{kF,uj) while the node would be detected in 
A{kp,uj). Thus, experimentalists should investigate the 
possibility of nodes at points in momentum space that do 



not belong to the normal state FS and they must keep 
in mind that some of the nodes may be determined by 
shadow bands with very small spectral weight. 

The most recent experimental results with the 
polarization dependence of the ARPES spectra for 
BaFei.85Coo.i5As2 provided the allowed contribution of 
each of the five 3d orbitals to the electron and hole Fermi 
surfacesi^ While discrepancies with proposed four and 
five orbital models are remarked in that publication, it is 
interesting to observe that the minimal two-orbital model 
addressed here does not contradict the ARPES findings 
if we disregard the additional /3 hole pocket FS that they 
present, which is reasonable because it has (1^2 _y2 char- 
acter, an orbital not included in the minimal model con- 
sidered here. In fact, along their T — M direction, which 
corresponds to our F — X, our first hole FS is purely d^z 
as it is their hole FS; our second hole FS, which arises 
upon folding our extended FS along X — Y is purely dyz , 
as it is their hole FS, and our electron FS is purely 
dyz as it is their 7 electron FS. Our second electron FS 
(obtained upon folding) has a purely d^z character, while 
their electron FS "f' /a' appears to be mostly dyz and d^y, 
but with some amounts of dxz as well. Along the diag- 
onal direction T — X, which corresponds to our F — M, 
the two hole Fermi surfaces are a symmetric admixture of 
dxz and dyz, exactly as in the two-orbital model.— Thus, 
these similarities between the experimental results and 
the band composition of the simple two-orbital model of- 
fer encouragement towards exploring whether the main 
physics of the pnictides can be captured with such a min- 
imum number of degrees of freedom. 

We have also mentioned in the text that the symmetry 
of the lattice and of the orbitals introduce constraints on 
the possible pairing operators. In general, a purely in- 
traband pairing, such as the proposed s± state,— '^'^'^ 
would occur only if the coupling of the electrons with the 
source of the attraction is identical for all orbitals. Con- 
sidering the different spatial orientations of the orbitals, 
it is not obvious that this should be the case, since, as 
discussed in the Introduction, phonons couple differently 
to electrons in the Pz and Px boron orbitals in the case 
of MgB2. Thus, if indeed the coupling results to be the 
same for all 3d orbitals we could use this fact to eluci- 
date the pairing mechanism; but, if this is not the case, 
we would expect some degree of interband pairing, at 
least in some regions of the Brillouin zone and, thus, it 
is important to study the experimental and theoretical 
consequences of such a possibility. 
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